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We show that the existence of the family of self-adjoint Lyapunov operators in- 
troduced in [J. Math. Phys. 51, 022104 (2010)] allows for the decomposition of the 
state of a quantum mechanical system into two parts: A past time asymptote, which 
' is asymptotic to the state of the system at t — > —oo and vanishes at t — oo, and a future 

time asymptote, which is asymptotic to the state of the system at t — >^ oo and vanishes 
f— I ' at t ^ oo. We demonstrate the usefulness of this decomposition for the description of 

I resonance phenomena by considering the resonance scattering of a particle off a square 

barrier potential. We show that the past time asymptote captures the behavior of the 
resonance. In particular, it exhibits the expected exponential decay law and spatial 
' probability distribution. 



I. INTRODUCTION 



pH , In standard non-rclativistic quantum mechanics time enters as a parameter, external to the quantum system being 
' investigated, and as such is not a dynamical variable associated in any way with the system's dynamics. However, 
^ I there are many cases, such as experiments measuring the time of arrival of particles at a detector, the decay time 
of an unstable quantum system etc., which call for a more dynamical point of view with respect to time. One 
approach to this problem consists of the construction of time operators through the use of covariant positive operator 
valued measures (POVMs) Since the time operators built via the use of POVMs are in general maximally 

^ ' symmetric, non-self-adjoint operators, the construction of these objects bypasses an old theorem of Pauli stating 
, that there does not exist a self-adjoint time operator T canonically conjugate to a Hamiltonian H whose spectrum 
' is semibounded, such that the pair T and H forms together an imprimitivity system However, in comparison to 
standard self-adjoint quantum observables these operators present some difficulties, e.g., in the relation between the 
algebra they generate and their spectral representations. 

The dynamical role of time within the framework of standard quantum mechanics has recently been considered 
in Q through the construction of a family of self-adjoint Lyapunov operators. Throughout this text we define 
a Lyapunov operator as a self-adjoint operator whose expectation value is monotonically decreasing in time. More 
precisely, let H be a Hilbert space corresponding to a given quantum mechanical system, let if be a self-adjoint 
Hamiltonian generating its evolution, and let \ip{t)) ~ U{t) ~ cxp{~iHt) \^) be the state of the system at time t. 
\ Define the trajectory ^f^, corresponding to an initial state \tp) £ H, as the set of states 

S : ^^4, {\m)}teR+ = iUit) mtm+ ■ (1) 
Then the definition of a Lyapunov operator is as follows Q : 

Definition 1: Let M be a bounded self-adjoint operator on %. Let 5*^, be a trajectory corresponding to an 

arbitrarily chosen initial state £ %. Denote by M(\E'^,) — {{ip\M\if) \ \ip) G the collection of all expectation 
values of M for normalized states in Then M is a forward Lyapunov operator if the mapping ti\j_.^ : R"*" i— >■ M('I'^) 
defined by 

TM,i.{t) = {m\M\m) (2) 

is monotonically decreasing in time. □ 

Remark 1: If in the definition above we require that tm,iP be monotonically increasing instead of monotoni- 
cally decreasing we also obtain a valid definition of a Lyapunov operator. The requirement that tm,-^ is monotonically 
decreasing is made purely for the sake of convenience. 
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It follows from this definition that a Lyapunov operator allows for the temporal ordering of states in any trajectory 
^",0 according to the ordering of the expectation values in M(^'^), thereby introducing temporal ordering into the 
Hilbcrt space % of any problem for which such an operator can be constructed. 

The construction of Lyapunov operators constitutes a somewhat conservative approach to the problem of time in 
quantum mechanics. In general, Lyapunov operators indicate the direction of flow of time in a system but do not serve 
as time operators and are not intended for answering questions of direct time measurements. Rather, the construction 
of Lyapunov operators, and in particular the family of Lyapunov operators introduced in Q, was carried out with 
a different goal in mind. The temporal ordering introduced into the Hilbcrt space is a fundamental property of a 
quantum mechanical problem admitting Lyapunov operators. A question then arises as to the possible implications 
this property may have on the description of the dynamics of the system, i.e. is there a way to 'inject' the direction of 
time ordering introduced by the Lyapunov operator into the description of the dynamics of a quantum system in such 
a way that certain processes in the evolution of the system obtain simple descriptions amenable to thorough analyses? 
The study of such consequences of the existence of a Lyapunov operator has been addressed to some extent in 0, Q . 
In the present paper we continue this line of investigation, and demonstrate the beginning of an affirmative answer to 
the question raised above. Specifically, we show that the existence of the family of Lyapunov operators introduced in 

Q leads to a certain description of the scattering process, which we shall term the transition decomposition, that 
is particularly useful for treating the evolution of scattering resonances. 



II. LYAPUNOV OPERATORS 



In this section we state some key results from js], which will provide the basis for the all that follows. The 
following theorem is an adaptation of results stated irL[5l, Q : 

Theorem 1: Let % he the Hilhert space representing some quantum mechanical system, let H be its Hamil- 
tonian, and let the spectrum of H be M"*", absolutely continuous, and uniformly degenerate. Then for any choice of 
generalized eigenhasis \E, A) of H, A being the degeneracy index, the operator 

1 foo poo 1 

where the summation over A may also stand for integration, is a Lyapunov operator in the sense of Definition 1, and, 
moreover, for any state ^"H 

lim (V- (t) |Af| (0) = , lim {i' {t)\M\i: {t)) = I . (4) 

t—¥oo t—^ — oc 



□ 

Remark 2: In this paper we shall consider only Hamiltonians for which the spectrum is R+, absolutely continuous, 
and uniformly degenerate. The results below apply also to the more general case of Hamiltonians for which the 
absolutely continuous spectrum is and is uniformly degenerate. In such cases the constructions below hold in the 
subspace Hac- 

Remark 3: Note that each complete set of generalized eigentates of H has associated with it a different Ly- 
paunov operator, and as such Theorem 1 represents a prescription for the construction of a family of self-adjoint 
Lyapunov operators. In particular, this still holds true for two sets that differ only by an energy dependent phase. 
When discussing general properties of the set of Lypaunov operators defined above, we shall often refer to M without 
specifying the underlying complete set of generalized eigenstates. 

Remark 4: The above family of self-adjoint Lyapunov operators has recently been generalized in Q- 

Remark 5: Throughout this paper we employ natural units, i.e. H = c = 1. 

The following corollary to Theorem 1 is crucial for what follows. 

Corollary 1: Let A A/^/^. Then for any state e U 

hm (^(t)|A|^(i)) =0, hm (V'(t)|A|^(0) = l. (5) 

□ 
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III. THE TRANSITION DECOMPOSITION 

111 this section we present a decomposition of the wave-function into two components - one that is asymptotic to 
the state of the system in the far past and vanishes in the far future, and one that is asymptotic to the state of the 
system in the far future and vanishes in the far past - and a corresponding decomposition for operators when working 
in the Heisenberg picture. The transition over time between the two asymptotic components suggests the possibihty 
that this decomposition yields an appropriate description of transient phenomena in the evolution of a quantum 
system. This expectation is realized when we apply the decomposition to scattering problems, where it will be seen 
to afford a highly useful description of the evolution of scattering resonances. This task is taken up in the next section. 

Let IV'(i)) = U{t) lip) € "H, t e M. At any moment in time t express ip{t) as the sum of two components as 
follows: 

m)) = \Mt)) + \^fit)), (6) 

where 

\Mt)) ■■= A m)) , \^f{t)) ■■= {I A) \m) ■ (7) 

The limits below are then an immediate consequence of Theorem 1 and Corollary 1: 

lim |||V;(i))-|V^6(t))|| =0, lim |||Vh(t))|| =0, (8) 
hm |||V/(t))|| =0, lim|||V'(t))-|V/W)|hO. (9) 

Eqs. (O, ([5]), and © imply that \ip{t)) can be decomposed into a sum of two components, \tpb{t)) and \tpf{t)), such 
that vanishes in the future time asymptote and is asymptotic to \ipit)) in the past time asymptote, and 

vanishes in the past time asymptote and is asymptotic to in the future time asymptote. We refer to li^bit)) 

as the backward asymptotic component and to {ipfit)) as the forward asymptotic component of \tlj{t)). Since the 
decomposition here is of the evolving state \'ip{t)), we call the decomposition in Eq. ^ the transition decomposition 
in the Schrddinger picture. In this decomposition the evolution of represented as a transition from its 

backward asymptotic component to its forward asymptotic component. 

The transition decomposition in the Schrodinger picture, Eq. ((S]), gives rise to a corresponding decomposition 
in the Heisenberg picture. Let X be a self-adjoint operator representing some physical observable and let X (t) = 
[t) XU (t) be its Heisenberg evolution. Consider the expectation value of X{t) for some arbitrary state E H 
and apply the decomposition in Eq. ([6]). We have 

(V' \x {t)\ ^) = (v- (t) \x\ i> {t)) = {i,b it) + i>f (i) \x\ {t) + V/ {t)) 
= (t) \x\ ^b W) + {{4'b W 1^1 V'/ it)) + {^f {t) \x\ v-fc (t))) + {^f (t) \x\ (t)) . (10) 

Define now the following decomposition of X , which will be seen below to correspond to the three different terms on 
the right-hand side of Eq. pH)) . 

X -.^Xb+Xtr+Xf , (11) 

with 

Xh-.^AXA, Xtr -.^ AX{I - A) + (I - A)XA, Xf:={I-A)X{I-A). (12) 
Indeed, if Xi,{t), Xtr{t), and Xf{t), respectively, are the Heisenberg evolutions oi Xt,, Xt,., and Xf, so that 

X{t) = Xb{t) + Xtr{t) + Xf{t) , (13) 

then 

\Xb {t)\^) = {t) ^ {t)) = (t) \AXA\ i' it)) = (^b {t) \X\ V^, (t)) , (14) 



(V- \Xtr (t)| V) = (V' (i) \Xtr\ ^ {t)) = (V' it) I AX (/ - A) + (/ - A) AA| {t)) 
= {^l,(t)\X\4,f{t)) + {^Pf{t)\X\i;b{t)) , 



(15) 
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\Xf {t)\ i,) = {i, (t) \Xf\ ^ it)) = {i, (t) \{I~A)X{I~ A)| V (t)) = (V-/ (t) \X\ 4,f (t)) . (16) 

Using the limits of |V'6(0) ^'^id |'(/;/(t)), Eqs. (jS]) and ([9]), we obtain that the components of X{t) satisfy the following 
limits 

\\m {X{t)~Xb{t))=Q, limX6(i) = 0, (17) 
lim Xtr(i) = , lim Xtr{t) = , (18) 
hm Xf{t)=Q, \un{X{t)-Xf{t))=0. (19) 

Eq. (|13p provides the transition decomposition in the Heisenberg picture. The observable X{t) decomposes into 
a sum of three components, Xb{t), Xtr(t), and Xf{t), such that Xb{t) vanishes in the future time asymptote and is 
asymptotic to X{t) in the past time asymptote, Xf{t) vanishes in the past time asymptote and is asymptotic to X{t) 
in the future time asymptote, and Xtr{t) is transient and vanishes both in the past and in the future time asymptotes. 
The evolution of X{t) is represented as a transition from Xi,{t) in the backward asymptote to Xf{t) in the forward 
asymptote. Accordingly, Xb{t) will be termed the backward asymptotic transition observable, Xf{t) will be termed 
the forward asymptotic transition observable, and Xtr(t) will be termed the transient observable. 

IV. APPLICATION TO SCATTERING PROBLEMS 

In this section we apply the transition decomposition to a quantum mechanical scattering problem. In particular, 
we find that this decomposition is especially useful for the description of scattering resonances. Specifically, we 
shall apply the transition decomposition to the evolution of approximate resonance states defined in the context of 
the formalism of the semigroup decomposition of resonance evolution Q. This formalism, developed in recent 
years out of efforts to adapt the Lax-Phillips scattering theory [l^] to the description of the evolution of scattering 



resonances in quantum mechanics |ll| . is a mathematical framework serving as a basis for a time-dependent theory 
of resonances in quantum mechanical scattering problems [l2j . We give here a short presentation of the simplest 
form of this framework that is sufficient for our purposes. 

Recall that Theorem 1 defines a large class of Lypaunov operators in the sense that for any complete set of 
generalized eigenstates of the Hamiltonian there corresponds a different Lyapunov operator. It follows that each 
choice of Lyapunov operator has associated with it its own transition decomposition. We would like to apply the 
transition decomposition to scattering problems. For such problems there are at least two distinguished energy 
representations for the Hamiltonian, i.e. the incoming and outgoing energy representations, and we must choose 
which of the two to employ. The incoming and outgoing energy representations are defined in terms of the incoming 
and outgoing solutions of the Lipmann-Schwinger equation. Let , A) and {E" , A) be the incoming and outgoing 
solutions of the Lipmann-Schwinger equation satisfying H\E^, A) = E\E^, A), where A stands for the degeneracy 
indices of the energy spectrum. (The degeneracy is assumed to be uniform over the energy spectrum.) The outgoing 
Lyapunov operator is defined as 

oo oo 

M+ = —YfdEfdE'\E',X) ^- -(E'-,X\, (20) 

-^0 

and gives rise to an outgoing transition decomposition in the Schrddinger or Heisenberg pictures. Defining A+ := 

1/2 

, the outgoing transition decomposition in the Schrodinger picture is given by 

\^it)) = \i'tit)) + H'fit))^ (21) 

while for an observable X defined on "H, in the outgoing transition decomposition in the Heisenberg picture, we have 

X = X++ X+ + X+ , (22) 

with 

X+ = A+XA+, X+^A+X{I~A+) + {I~A+)XA+, X+ ^ {I - A+)X{I ~ A+) . (23) 
To proceed we assume the following: 
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(i) Let Hhe a Hilbert space corresponding to a given quantum mechanical scattering problem. A self-adjoint 'free' 
unperturbed Hamiltonian Hq and a self-adjoint perturbed Hamiltonian H are defined on Ji and form a complete 
scattering system, i.e., we assume that the M0ller wave operators il^{Ho, H) exist and are complete. 

(ii) We assume that the (absolutely) continuous spectrum of H is uniformly degenerate. To simplify matters we 
assume that this degeneracy is one. 

(iii) The S'-matrix in the the energy representation, denoted by S{E), is the boundary value of a function S{z) 
analytic in some strip above the positive real axis and having an analytic continuation across the cut on the 
positive real energy axis into some simply-connected region E below the real axis in which it has a single, simple 
resonance pole at the point z = /i with Im/i < 0. 

It is shown in [1,0] that there exists a dense set S C H and a well defined state \ip'^^^) € T-L such that for any states 
\ip) G E and \ip) G % the above assumptions lead to a decomposition, induced by the pole of the ^-matrix at z = /x, 
of matrix elements of the evolution U {t) generated by H, having the form 

{^\uim=Bi^A',^^■,t) + a{^,^l){r^\^)e~'''*^ ^>o. (24) 

The second term of the decomposition on the right-hand side of Eq. (|24p exhibits the typical decay behavior of a 
resonance. This term is the semigroup term, or resonance term. The term B(ip,ip, fi]t) in Eq. p4p is the so-called 
background term. The state \tjj'^^^) appearing in the second term in Eq. (j24p is called the approximate resonance state. 
Note that if the state \tjj) on the left-hand side of Eq. ((24|) is chosen to be orthogonal to lip"^^^) then the resonance 
term in that equation vanishes. The reference to \tp'^^^) as an approximate resonance state stems from the fact that it 
can be shown that there is no choice of \ip) and \ip) in the matrix clement {f \ U(t)ijj) for which the background term 
B{if,tp,iJ.;t) disappears [1]. Indeed, the Schrodinger evolution of a closed system docs not allow for an exponential 
decay law for the survival probability [l^ and deviations, such as the Zeno effect for short times arc inevitable [l^ . 
An explicit expression for the approximate resonance state ]?/'^^^) is given by [1, H, HB] 

This expression is obtained under the assumption that there is only a single resonance pole of the 5'-matrix below 
the positive real axis in the region S. In the case that E contains multiple resonance poles of the S'-matrix Eq. ([25)1 
provides only a zeroth order approximate resonance state [^. Throughout the rest of the paper we only consider 
approximate resonance states given by Eq. (|25p . This restriction pertains also to the example worked out below for 
which the S'-matrix possesses multiple resonance poles. Hence, all states calculated there are of zeroth order. The 
restriction to zeroth order approximate resonance states is made for the sake of simplicity and clarity of exposition 
and it should be emphasized that there is no a priori difficulty in working with higher order approximate resonance 
states. 

Applying the outgoing transition decomposition in Eq. ([22]) to the expectation value of an observable X(t) in the 
state \i>'^PP) 1| IV'^^) H"^ l^^^), we get 



(v;nx(i)i<^o = (v^r + {r,'"'\xut)\^r) + (v;rix+(t)i^7^') . (26) 

We shall use this decomposition to represent the evolution of a particular resonance in a simple one-dimensional 
scattering problem. The model we consider is the scattering along the half-line K"*" off a square barrier potential. 
Thus, we consider a free Hamiltonian Hq = —^d^ acting on L^(M+) (where Hq is taken to be the self-adjoint 
extension in L'^{R+) of -^d^ from its original domain of definition V{dl) = {0(a;) ] 0(a;) e Wi{M.+ ), (j){0) = 0}) 
and a full Hamiltonian H = Ho + V with V a multiplicative operator {Vip){x) = V{x)ip{x) such that 

{0 , < a; < a , 
Vo , a<x<b, (27) 
0, b <x, 

where < a < 6 and Vq > 0. In this case there are no bound state solutions of the eigenvalue problem for H and the 
(absolutely) continuous spectrum of H is K+. In order to find the scattering states and calculate the S'-matrix one 
solves the eigenvalue problem for the continuous spectrum generalized eigcnfunctions ipEix) of H 

(^-^dl + Vix)^i;E{x)=EijE{x), EeR+. (28) 
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ir^^i^f [m] 




|(^;,f (E)p [m-'] 
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E[m] 



Figure 1: Spatial density of the third resonance fj,3- 



Figure 2: Energy density of the third resonance //a- 



Imposing appropriate boundary conditions wc find that 



iPe{x) 



ai (fc) sin kx 



h/32(fc) 



< X < a, 
a < X < b , 
b < X , 



(29) 



where k = y/2mE and k' = y/2m{E - Vq) for E > Vq > or k' = iV2mVb - 5^ for Vb > £^ > 0. The coefficients in 
Eq. (gH) are given by ^ 



a2(k) = -e ° I sinfca + cosfca ) ai(fc) 
2 \ ik' J 



sin ka — — cos ka ai (k) 



asik) = ie-"^-^ ( fl + ^)e*^-'(^-'^)(sinfca+ A cosfca) 
4 \\ fc/ ifc' 



(30) 



^' ')e-*'='(''-'^) (sin fc-a - cosfca j ) ai(fc) 



/33(fc)-a3(fc) > 

where ai(fc) is to be determined by normalization conditions. 

Given the full set of solutions {'ipEix)}E£R+ for the continuous energy spectrum one can find the sets }_e6r+ of 
solutions of the Lippmann-Schwinger equation corresponding to incoming and outgoing asymptotic conditions. We 
have 



V'+(x) {x\E+) = 



2a* (fc) 



V'b(x) 



(x\E- 



1 



2^3 (fc) 



■)pEix) 



(31) 



where ^/'^(a;) and '4)^{x) are, respectively, the incoming and outgoing Lippmann-Schwinger solutions. The normaliza- 
tion conditions for the Lippmann-Schwinger states in Eq. (j3ip give us ai (fc) = (27rfc)^^/^. In the energy representation 
the iS-matrix is given by 



S{E) 



«3(fc) 

a^(fc) 



(32) 



The above expression for the S'-matrix leads to the calculation of the scattering resonances of the problem. For a 



E, 



iT^j2 with E^^ > 



resonance point z = fxj in the lower half-plane below the positive real axis we set jij 
the resonance energy and F^^ > the resonance width. 

For barrier parameters a — 2m~^, b = 3m~^, and Vq = 5m the three lowest energy resonance poles arc given 
by fii ~ 0.9106 m- i 0.0012 m, fi2 ~ 3.51 19 to- i 0.0282 m and ^3 7. 1168 to- 1 0.4462 to. We shall focus on the 
third resonance pole ^3. Utilizing the outgoing Lippmann-Schwinger eigenfunctions tjj^{x), given by Eqs. (|29ll3ip . 
the spatial wave function of the approximate resonance state ijj'j^P^{x), and its energy density, can be calculated 
numerically. The probability density I'fp'^^^ {x)\'^ is shown in Fig. 1 (see Q), while the energy density I'lp'/i^^ {E)\'^ is 
shown in Fig. 2. 
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Figure 3: Monotonic decrease of the expectation value of in the state \'ip1^{t)). Note the transition to exponential decay 
for t > 0. 



Now apply the outgoing transition decomposition of Eq. (|26p to the evolution of the approximate resonance state 
\i(j'^PP{t)) = U{t)\ipi;P/). According to Eq. ^ the spatial probability distribution of lifj^f/ {x, t) p decomposes into 
three components 

i<r or = \4'ZT or + 2Re((o^ o)*<:y {x, o) + {x, or • (33) 

The right-hand side of Ec^. (p3|) is the outgoing transition decomposition of the position probability density of 'ip'^^^it)- 
The first term on the right hand side of Eq. p3p is the backward asymptotic component, the second term is the 
transient component, and the third term is the forward asymptotic component. 

Fig. 4 presents the results of the application of the outgoing transition decomposition to the evolution of the 
approximate resonance state IV'^g''). Each row represents a 'snapshot' corresponding to a particular time t. The left 
graph in each row shows \i^'^P{x, i)p at time t. The middle graph shows the contribution of its backward asymptotic 
component |'0^3^h^(a;, OP *t time t, and the right graph in each row shows the sum of transient and forward asymptotic 

components 2Re((i/j°||'^^^(x, t))*'ip'ZJ''f^ [x, t)) + lV')^3^/^(2;, It is clearly seen from this sequence of snapshots that 
the formation phase of the resonance, starting at the negative time asymptote t — >■ —oo and following through the 
scattering process up to t = 0, is captured by the backward asymptotic component (the middle column of graphs), 
while through this whole time interval the contribution of the transient and forward asymptotic components is small. 
In the decay phase of the resonance, commencing at t > 0, the middle column of graphs essentially shows the spatial 
probability density profile of a resonance state, multiplied by an exponentially decaying factor exp(— F^gt), which 
gives the decay of the resonance over time. As the resonance state decays the probability is transferred to the forward 
asymptotic term, a process captured in the right most column of graphs, and is eventually carried to spatial infinity 
as the scattering process evolves further towards the forward time asymptote i — >■ oo. Observe also that the graphs 
on the right column in Fig. 4 obtain also negative values. This is due to the contribution of the transient term which 
is not necessarily positive. 

Further understanding of the behavior of the probability density \'4''Z^'b^ {^^t)]'^ may be gained by integrating it 
over X. In this way we get 

oo oo 

j dx t)\' ^ J dx{rzT it) I x){x I v^r.^ (0) 



oo 

= J dx{A+r,7 it) I x){x I A+r,7 it)) = (O v^^f (t)) , (34) 



which is just the expectation value of the Lyapunov operator in the state (t)) ■ A plot of this behavior 

is presented in Fig. 3. The decay of the expectation value of A/+ matches the resonance decay factor exp(— F^gf) 
mentioned above. 
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Figure 4: Transition decomposition of ■ip'j^^^(t). Each row of graphs represents a time frame at time t indicated on 
the left row. The left row depicts {x, t)\'^ , the middle row depicts lip""^^'^^ {x, t)\-^ , while the right row depicts 

\i^;i';-+(x,tr. 
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